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This Letter describes the current most precise measurement of the W boson pair production 
cross section and most sensitive test of anomalous WWy and WW Z couplings in pp collisions at a 
center-of-mass energy of 1.96 TeV. The WW candidates are reconstructed from decays containing 
two charged leptons and two neutrinos, where the charged leptons are either electrons or muons. 
Using data collected by the CDF II detector from 3.6 fb _1 of integrated luminosity, a total of 654 
candidate events are observed with an expected background contribution of 320 ± 47 events. The 
measured total cross section is o~{pp — > W + W~ + X) = 12.1 ± 0.9 (stat) (syst) pb, which is 
in good agreement with the standard model prediction. The same data sample is used to place 
constraints on anomalous WWj and WW Z couplings. 

PACS numbers: 14.70.Fm, 13.38.Be, 13.85.Qk 
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The measurement of W boson pair production is an 
important test of the standard model (SM) of particle 
physics. This process is also an essential background to 
understand for Higgs boson searches at particle collid- 
ers. Next-to-leading order (NLO) calculations of W + W~ 
production in pp collisions at \fs — 1.96 TeV predict 
a cross section of a NLO {pp -> W+W~) = 11.7 ± 0.7 
pb 0, 0] • The presence of anomalous WWj and WWZ 
triple gauge boson couplings (TGC) [j| could be indica- 
tions of new physics at a higher mass scale, and would 
lead to rates for W + W~ production or kinematic distri- 
butions that differ from those predicted by the SM. 

This Letter reports a measurement of the W + W~ pro- 
duction cross section and limits on anomalous TGCs us- 
ing a final state consisting of two oppositely charged lep- 
tons and two neutrinos in pp collision data collected by 
the CDF II detector from 3.6 fb _1 of integrated lumi- 
nosity. First evidence for W boson pair production was 
reported by CDF using Tevatron Run I data This 
process was later measured with greater significance by 
CDF and DO using 184 pb" 1 and 224-252 pb^ 1 respec- 
tively of integrated luminosity from Run II of the Teva- 
tron [10. Recently, the DO collaboration measured the 
W + W~ cross section with a precision of 20% using 1.0 
fb _1 of integrated luminosity from Run II [7j. Limits on 
anomalous TGCs have previously been reported by LEP 
experiments and the CDF and DO collaborations p], @| • 

The cross section measurement uses a matrix element 
method in which the probability for each event to have 
been produced by several relevant SM processes is cal- 
culated. A likelihood ratio (LR) is formed from these 
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probabilities. The predicted shapes and normalizations 
of the signal and background LR distributions are used 
to extract the SM W + W~ production cross section via a 
maximum-likelihood fit to the LR distribution observed 
in data. In general, the presence of anomalous TGCs will 
increase the number of events containing leptons with 
very high values of momentum transverse to the beam 
direction, as compared to SM expectations. Limits on 
anomalous TGCs are determined from the shape and nor- 
malization of the transverse momentum spectrum con- 
structed from the highest transverse momentum lepton 
in the event, referred to as the leading lepton. The results 
are reported in the HISZ scheme, where three parame- 
ters, \z, gf , and K ry , ell"© used to describe all dimension- 
six operators which are Lorentz and SU(2)j,® U{\)y 
invariant and conserve C and P separately [9j. In the 
SM, \z — and gf = k 7 = 1. In this Letter, Agf and 
Ak 7 are used to denote the deviation of the gf and k 7 
parameters from their SM values. The non-SM values 
of the parameters Xz, g± , and k 7 are functions of the 
invariant mass of the W + W~ system, \/~§. These results 
probe a larger range of values of yfi and thus may only 
be qualitatively compared to the results from LEP, which 
were below \/§ = 209 GeV. For hadron collisions, a dipole 
form factor for an arbitrary coupling a(s) = ^ 1+ "° A p 

is introduced to turn off the coupling at large y/J and 
avoid a violation of unitarity. The form factor scale A is 
the scale of new physics. 

In the CDF II detector [1(3], a particle's direction is 
characterized by the azimuthal angle <j> and the pseudora- 
pidity r\ = — ln[tan(0/2)], where 9 is the polar angle mea- 
sured from the proton beam direction. The transverse 
energy Et is defined as Es'm6, where E is the energy 
in the calorimeter towers associated with a cluster of en- 
ergy deposition. Transverse momentum, px, is the track 
momentum component transverse to the beam line. The 
magnitude of the pt for an electron is scaled according 
to the energy measured in the calorimeter. The missing 
transverse energy vector, $ T , is defined as — J2i ™t-> 
where h T is the unit vector in the transverse plane point- 
ing from the interaction point to the energy deposition in 
calorimeter tower i. This is corrected for the px of muons, 
which do not deposit all of their energy in the calorimeter, 
and tracks which point to uninstrumented regions in the 
calorimeter. The scalar $ T is defined as \$ T \. Strongly 
interacting partons produced in the pp collision undergo 
fragmentation that results in highly collimated jets of 
hadronic particles. Jet candidates are reconstructed us- 
ing the calorimeter and are required to have Et > 15 
GeV and \rj\ < 2.5. Isolated lepton candidates are ac- 
cepted out to an 1 77 1 of 2.0 for electron candidates and |ry| 
of 1.0 for muon candidates. 

The experimental signature for the decay W + W~ — > 
£ + v£~D is two reconstructed leptons with opposite charge 
and l£ T from the neutrinos which escape undetected. In 
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this Letter, £ refers to an electron or muon. Additional 
signal acceptance (~12%) is obtained from cases where 
one or both W bosons decay to a tau lepton which sub- 
sequently decays to an electron or muon. There are sev- 
eral SM processes which result in a similar final state to 
W + W~ and are therefore backgrounds in this measure- 
ment. These are other diboson production (WZ, ZZ) 
and top-quark pair production (tt). It is also possible to 
observe apparent ]ft T arising from the mismeasurement of 
lepton energy, lepton momentum, or the hadronic part of 
the final state. Drell-Yan (Z/j* — > £ + £~) events have no 
neutrinos in the final state, but due to large production 
rates enter the W + W~ candidate sample via mismea- 
surements. A third source of background is events in 
which a final-state particle is misidentified. These are 
IU+jets and W*y production, where the W boson decays 
lcptonically and a jet is reconstructed as a lepton can- 
didate or the 7 converts in the detector material and is 
reconstructed as an electron. 

Events containing two oppositely charged lepton candi- 
dates are selected from the data sample. The online event 
triggering and selection of lepton candidates are identical 
to those used in the search for SM Higgs bosons decay- 
ing to two W bosons at CDF 1 lj . The leading-lepton px 
is required to be above 20 GeV/c to satisfy the trigger 
requirements, while the second lepton is allowed to have 
a pt as low as 10 GeV/c. The requirement is also made 
that events contain no jet candidates, which significantly 
reduces tt background. A variant of $ T used in selecting 
candidate events is defined as $ T rel = ]fi T sin A(f>($! T , £) 
when A<j>($ T ,l) < where A<j)(]ft T ,£) is the azimuthal 
separation between the and the momentum vector of 
the nearest lepton candidate. If Acf)(]p T ,£) > then 
$t rei — $t- The lfi T rel variable is designed to reject 
events where the apparent $ T arises from the mismea- 
surement of lepton energy or momentum, and is required 
to be above 25 GeV to reduce the otherwise large Drell- 
Yan contamination. This requirement is lowered to 15 
GeV for electron-muon events where contributions from 
Drell-Yan are inherently smaller. The Wj and heavy- 
flavor (J/ip,Y) backgrounds are reduced by requiring 
that the invariant mass of the lepton pair be greater than 
16 GeV/c 2 . 

With the exception of the W + 1-jet background, 
the acceptance and kinematic properties of the signal 
and background processes are determined by simulation. 
Events from W + W~ are simulated at NLO using the 
mc@nlo generator @. The tt, WZ, ZZ, and Drell- 
Yan backgrounds are simulated with the PYTHIA gener- 
ator [13]. The Wj background is determined using the 
generator described in Ref. The response of the 

CDF II detector is modelled with a GEANT-3-based sim- 



TABLE I: Expected number of signal (W + W~) and back- 
ground events along with the total number of expected and 
observed events in the data. Uncertainties include all system- 
atic uncertainties described in the text. 



Process 


Events 


Z/7* (Drell-Yan) 


79.8 


± 


18.4 


WZ 


13.8 


± 


1.9 


W-y 


91.7 


± 


24.8 


W + 1-jet 


112.7 


± 


31.2 


ZZ 


20.7 


± 


2.8 


tt 


1.3 


± 


0.2 


Total Background 


320.0 


± 


46.8 


w + w~ 


317.6 


± 


43.8 


Total Expected 


637.6 


± 


73.0 


Data 


654 





ulation [14[ . The expected yields for each process are nor- 
malized to the cross sections calculated at partial next- 
to next-to-leading order (ft 15]), NLO (W+W~ 0, 0, 



WZ and ZZ [l|), or leading-order with estimated higher- 
order corrections (Wj [l3[ and Drell-Yan [3]). Effi- 
ciency corrections for the simulated detector response to 
lepton candidates are determined using samples of ob- 
served Z — > £ + £~ events. The W + 1-jet background is 
calculated using the probability, measured in indepen- 
dent jet-triggered data samples, that a hadronic jet will 
be reconstructed as a lepton candidate. These probabil- 
ities are applied to the jet in the W + 1-jet data sample 
to estimate the number of such events which will pass 
the full lepton identification and signal selection crite- 
ria. The expected signal and background contributions 
are given in Table U along with the observed number of 
events. 

The dominant systematic uncertainties on the esti- 
mated contributions come from the luminosity measure- 
ment (6%) [13] and the simulated acceptances of the sig- 
nal and background processes. The acceptance uncer- 
tainty due to the parton distribution function modeling 
ranges from 1.9% to 4.1% for the different processes. A 
10% uncertainty is assigned to all simulated processes 
for the kinematic differences between leading-order and 
higher-order calculations, based on the difference in ac- 
ceptance of W + W~ events simulated at leading-order 
and NLO using the pythia and mc@nlo generators re- 
spectively. This uncertainty is reduced to 5% for the 
W + W~ signal, which is simulated at NLO. The cross 
section uncertainties are 6% on diboson production, 10% 
on ft and W'y production, and 5% on Drell-Yan produc- 
tion. A 21% uncertainty is included for the Drell-Yan 
background to account for the mismodeling of Tp T and 
jet production rates. Systematic uncertainties of 20% 
and 27% are assigned to the W'y and W + 1-jet back- 
ground estimates, respectively, due to uncertainties in 
the modeling of the photon conversions and misidentifi- 
cation of a jet as a lepton. Uncertainties on the modeling 
of jets accounts for 2% to 4% and lepton identification 
and trigger efficiencies range from 1% to 7%. 

For each event passing the signal selection criteria, four 
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matrix-element-based event probabilities are calculated 
corresponding to the production and decay processes 
W+W- -► e+v£~v, ZZ -► t+l-vv, W + 1-jet -> 
1-jet, and W7 — * £^ + 7. In the latter two processes, 
the jet or 7 is assumed to have been reconstructed as a 
charged lepton candidate. The event probability for a 
process X is given by 



Px(x) 



1 



da(y) 



e(y) G{x, y) dy 



(1) 



where x represents the observed lepton momenta and 
Ifirp vectors, G(x, y) is a transfer function representing 
the detector resolution, and e(y) is an efficiency function 
parametrized by r\ which quantifies the probability for a 
particle to be reconstructed as a lepton. The differen- 
tial cross section ^tf is calculated using leading-order 
matrix elements from the mcfm program [1] and inte- 
grated over all possible true values of the final state par- 
ticle 4-vectors y. The normalization factor (a) is deter- 
mined from the leading-order cross section and detector 
acceptance for each process. These event probabilities 
are combined into a likelihood ratio 



LR 



WW 



WW 



P\ 



WW 



J2j kjPj 



(2) 



where j = {ZZ, W + 1-jet, W7} and kj is the relative 
fraction of the expected number of events for the j-th pro- 
cess such that J^j kj = 1- The templates of the LRy/w 
distribution are created for signal and each background 
process given in Table [U 

A binned maximum likelihood is used to extract the 
W + W" production cross section from the shape and 
normalization of the LRyyw templates. The likelihood 
is formed from the Poisson probabilities of observing rij 
events in the z-th bin when ^ are expected. Variations 
corresponding to the systematic uncertainties described 
previously are included as normalization parameters for 
signal and background, constrained by Gaussian terms. 
The likelihood is given by 



i', 



where 



E 



{N E XP) 



(3) 



(4) 



f% is the fractional uncertainty for the process k due to 
the systematic c, and S c is a floating parameter asso- 
ciated with the systematic uncertainty c. The correla- 
tions of systematic uncertainties between processes are 
accounted for in the definition of The expected num- 
ber of events from process k in the i-th bin is given 
by (Vfc p ) t . The parameter au is an overall normal- 
ization parameter for process k and is fixed to unity 




~2 SM W + W" 
P^j Background 

-\- Data 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Matrix Element Likelihood Ratio (LR ) 

WW 



FIG. 1: The LRww distributions for the signal (W + W~) 
and background processes after a maximum likelihood fit to 
the data. 



for all processes other than W + W~ , for which it is 
freely floating. The likelihood is maximized with re- 
spect to the systematic parameters S c and aww us- 
ing the minuit program [l8j]. The W + W~ cross sec- 
tion is then given by the fitted value of aww multiplied 
by a NLO (pp — > W + W~). This method gives a mea- 
sured value for the W + W~ production cross section of 
a(pp -> W+W- +X) = 12.1 ± 0.9 (stat) t\% (syst) pb. 
The fit to the data of the signal and sum of the individ- 
ually fitted background templates is shown in Fig. [TJ 

The likelihood of the observed leading-lepton pt dis- 
tribution is used to set limits on anomalous TGC values. 
The robustness of the leading-lepton pt distribution has 
been verified using the same lepton selection in several 
non-overlapping final state kinematic regions. The re- 
sponse of the detector to events with different coupling 
constants is simulated for six points in the parameter 
space near the existing limits [8[. The efficiency multi- 
plied by acceptance as a function of the leading-lepton px 
is taken to be the average of the values measured in these 
samples. The uncertainty is taken to be the maximum 
variation among these samples and ranges from 7% at low 
Pt to 50% at high pr- This pr-dependent efficiency is ap- 
plied to the NLO generator-level distributions produced 
by the mcfm program 1] to predict the leading-lepton 
Pt spectrum for the coupling values considered, as shown 
in Fig. H 

Each of the likelihoods C(Xz), £(Agf ), and £(Ak 7 ) 
are computed as the product over all bins in the leading- 
lepton pt distribution of the Poisson probability of each 
bin given the model, and 95% confidence levels are set 
where (-21n£) - (-21n£ min ) = (1.96) 2 . The sys- 
tematic uncertainties include all those described for the 
W + W~ cross section and the additional pr-dependent 
uncertainty on the efficiency described previously. Sys- 
tematic uncertainties are implemented by simultaneously 
applying all variations which reduce the sensitivity. The 
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Lepton p T [GeV c" 1 ] 
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FIG. 2: Leading-lepton pr distribution for data compared to 
the SM expectation. Also shown is how the expectation would 
be modified by anomalous couplings near the observed limits. 

TABLE II: Expected and observed limits on anomalous 
TGCs. For each coupling limit set, the two other couplings 
are fixed at their SM values. Values of the couplings outside of 
the given observed range are excluded at the 95% confidence 
level. 





A (TeV) 


Xz Ag? 




Expected 
Observed 


1.5 
1.5 


(-0.05,0.07) (-0.09,0.17) 
(-0.16,0.16) (-0.24,0.34) 


(-0.23,0.31) 
(-0.63,0.72) 


Expected 
Observed 


2.0 
2.0 


(-0.05,0.06) (-0.08,0.15) (-0.20,0.27) 
(-0.14,0.15) (-0.22,0.30) (-0.57,0.65) 



observed 95% confidence limits, shown in Table [Til are 
weaker than expected. The probability of observing these 
limits in the presence of only standard model W + W~ 
production ranges from 7.1% to 7.6% depending on the 
coupling constants {\z,9i >k 7 ). 

In summary, the W + W~ production cross section has 
been measured in pp collisions at ^/s = 1.96 TeV from 
reconstructed events in the dilepton final state using a 
likelihood ratio formed from matrix-element-based event 
probabilities. This result is the most precise measure- 
ment at this energy with an overall uncertainty of less 
than 15%. The same event sample is also used to per- 
form the most sensitive probe to date at this energy of 
anomalous WW Z and WW"f couplings. The leading- 
lepton pt distribution of the sample is found to be in 
moderate agreement with the SM expectation and used 
to place limits on anomalous triple gauge couplings. 
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